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Abstract 

The convergence of the barrier height for the O -t- H 2 — » OH + H reaction has 
been studied as a function of the size of the active space and basis set completeness. 
The barrier height is rapidly convergent with respect to expansion of the active 
space. Addition of 2p— >2p' correlation terms to the active space lowers the barrier 
to the O -f H 2 reaction by about 2.0 kcal/mole, but addition of 3d and other terms 
has little additional effect. Multireference singles and doubles contracted Cl plus 
Davidson’s correction calculations using a [5s5p3d2flg/4s3p2dlf] basis set with a 
5o2tt active space lead to a barrier height of 12.7 kcal/mole. Including an estimate 
of the Cl contraction error and basis set superposition error leads to 12.4 kcal/mole 
as the best estimate of the barrier height. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 


1 



I. Introduction 

The national aerospace plane(NASP) will be an air breathing hypersonic flight 
vehicle capable of achieving velocities as high as Mach 25.0. It will utilize supersonic 
combustion ramjet( scramjet) engines and burn hydrogen for fuel. The high flow 
velocities introduce design difficulties, since the flame propagation speed may not 
be as large as the flow velocity i.e. the flame may blow out. This has resulted in 
considerable interest in the kinetics of reactions involving H, O, and N. Recently 
a research program has been started to compute from first principles the rates 
of selected reactions important in SCRAM jets. The computation of the rate of 
these reactions first entails the computation of the potential energy surfaces for the 
reactions of interest. 

In order to assess the accuracy of potential energy surfaces involving oxygen, a 
systematic study was carried out of the convergence of the barrier height for the 
O + H 2 — * OH -f H reaction as a function of both the size of the active space in 
the CASSCF calculation and the size of the basis set. This reaction was selected 
for study since a number of dynamics studies[l-4j have been carried out based on 
the POL-CI surface of Walch et al. [5] . The agreement between the computed and 
the experimental thermal rate constants establish that the barrier height of 12.5 
kcal/mole from the POL-CI surface is accurate to better than 1.0 kcal/mole. We 
also note here that a similar study by Howard et al.[6] using a first order Cl iterative 
natural orbital procedure led to a barrier height of 14.0 kcal/mole. 

Although the POL-CI barrier height is believed to be quite accurate, the com- 
puted endoergicity(AE ri ) on the other hand is -0.1 kcal/mole compared to an ex- 
perimental value of -(-2.9 kcal/mole. Thus, one has to assume that the good result 
for the computed barrier height is somewhat fortuitous. As discussed in Section 
III., attempts to improve upon the POL-CI result by using multireference singles 



and doubles Cl methods (MRSD Cl) lead to a barrier of 16.0 kcal/m<?le using the 
same three reference configurations and a comparable basis set to that used in the 
POL-CI study. Similarly, for the F + H 2 reaction a first order Cl calculation by 
Bender and coworkers[7] gave a barrier height of 1.66 kcal/mole, but more recent 
calculations using MRSD Cl and perturbation theory methods have lead to barriers 
greater than 3.2 kcal/mole even when very large basis sets are used [8]. 

In the present paper we examine the source of the differences between the POL-CI 
and MRSD-CI calculations. For both of these reactions, ionic terms are expected 
to be important in the saddle point and products regions but not in the reactants 
region. Thus, the problem of describing the barrier height for these reactions is 
closely related to the problem of dscribing the electron affinity of the O and F atoms. 
Botch and Dunning[9] found that 2p-*2p' terms were important in describing the 
O and F atom electron affinities. For example, with the largest basis set used in 
the present studies, the computed electron affinity for O is 1.27 eV for a MRSD-CI 
calculation based on a CASSCF including 2p — * 2p' correlation. This number may 
be compared to 1.46 eV for experiment [lOj and 1.04 eV for SCF plus singles and 
doubles Cl with a very large basis set[ll|. Thus, it is seen that 2p — *• 2p' correlation 
is very important in describing the electron affinity of O atom, and these terms were 
added to the active space in these calculations. 

Section II describes the basis sets and the computational approach. Section III 
presents the results of calculations using an expanded active space and extended 
basis sets, while the conclusions are presented in Section IV. 

II. Computational Details. 

Table I gives the O basis sets which were used in this study. The O sp basis set is 
a (5s3p] segmented contraction of the van Duijneveldt (lls6p) primitive set[l2j. The 
inner nine primitive functions were contracted (63) based on the Is orbital, while 
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the outer five s functions were contracted (311) based on the 2s orbital. (Note that 
this contraction uses three functions twice.) In order to describe O - character, the 
basis set was augmented by a single set of 2p functions with the exponent selected 
using an even tempered criterion leading to a [5s4p] sp basis set. The polarization 
function basis sets were developed from the gaussian fits to Slater type orbitals 
given by Stewart [ 13]. In each case the effective Slater exponent was optimized at 
the Cl level using a two term GTO fit[l3j to a Slater. The exponent for the 03d 
function was optimized at the Cl level for the ground state of the O atom, while the 
exponents for the 04f and 05g were optimized at the Cl level for the OH molecule. 

The H basis set is given in Table II. The hydrogen s basis set used the van 
Duijneveldt 6s set[l2j contracted to [3s] plus one diffuse function selected in an 
even tempered fashion. The [2p] basis set was selected as a four term fit|l3] to a 2p 
STO with exponent of 1.0. This primitive basis was contracted (211) and the most 
diffuse function was discarded leading to a [4s2p] basis set, which was used for the 
smallest basis set for O + H 2 . The scale factor of the H p function was subsequently 
optimized at the SDCI level for H 2 using a [4slp] basis set. This lead to an effective 
2p Slater exponent of 1.84. The [4pj basis was generated from this Slater exponent 
as a 211 contraction of a 4 term ST0[l3j plus one diffuse 2p function selected by 
an even tempered criterion. A [3p' set was also used which consisted of the [4p] set 
minus the most diffuse exponent. The 3d and 4f exponents used here were optimized 
at the Cl level for OH. From these s, p, d, and f basis functions five different basis 
sets were generated: [4s2p], [4s3pldj, [4s4pld], [4s4p2dj, and [4s4p2dlf]. 

During the course of this work, Almlof and Taylor implemented a new ver- 
sion of the molecule integral program which permits general contractions. These 
authors also developed basis sets for use with this program in which the occu- 
pied and virtual orbitals are generally contracted based on the natural orbitals 
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from an atomic Cl calculation[l4]. These basis sets are optimal for describ- 
ing the atomic correlation and thus have very small basis set superposition er- 
rors. In addition the basis sets are found to be sufficiently flexible to be used 
in molecular calculations at both the SCF and Cl level. In the present studies 
a (I3s9p6d4f2g/8s6p4d3f)/[5s5p3d2flg/4s3p2dlf] basis set[l4] was used. Here the 
fifth 02p function was a diffuse function to describe 0~ character. This last func- 
tion was selected using an even tempered criterion. This basis set is referred to as 
the atomic natural orbital (ANO) basis set. 

Most of the calculations were carried out on the CRAY XMP-48 with the 
MOLECULE!l5j-S\VEDENjl6j system of programs. Some of the basis set tests were 
carried out on the Cyber-205 using the Karlsruhe adaptation [17] of the COLUM- 
BUS [18] codes. These calculations used the Coupled Pair Functional (CPF) method 
of Ahlrichs[l9j. 

Since the Cl calculations reported here use a limited set of reference configura- 
tions, the Cl energy is not invariant to a unitary transformation of the CASSSCF 
orbitals among themselves. The CASSCF active orbitals are uniquely defined by 
a natural orbital transformation of the converged active orbitals, followed by a 
CASSCF Cl in the natural orbital basis to define the important reference configu- 
rations. 

Most of the Cl calculations presented here were uncontracted; however, some cal- 
culations were carried out with the externally contracted Cl of Siegbahn[20’ r Cal- 
culations which use the contracted Cl are denoted by CCI while the uncontracted 
calculations are denoted by Cl. In all cases the multireference analog of Davidson’s 
correction [21 ] was computed and results are reported both with and without David- 
son’s correction. The correction used in the CCI is A E (1 - Cq)/C£ which differs 
slightly from the original Davidson’s correction, used in the uncontracted Cl, which 
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is AE (1 — Cp), where AE is the Cl energy minus the reference energy and C 2 is 
the square of the coefficient of the reference configuration or the sum of the squares 
of the coefficients of the reference configurations in the multireference case. The 
results where Davidson’s correction has been applied are denoted by Cl + Q or CCI 
-I- Q for the uncontracted and contracted Cl, respectively. All of the calculations 
were performed for collinear geometries using C 2 v symmetry. 

III. Discussion. 

For collinear geometries, the SCF configuration of the reactants O + H 2 on the 
3 H surface is: 

lo 2 2a 2 3a 2 4a J l7r 2 l7r* (1) 

Here la and 2a are Ols and 02s like orbitals. The 3a orbital is the H 2 bonding 
orbital (1 a g ) and the 4a and l7r orbitals are the 02p orbitals. As described else- 
where[5], the 3a and 4 a orbitals map into the OHo bonding orbital and a His orbital 
for the products. In order to describe the spin recoupling process it is necessary to 
also include the configurations. 

Ia 2 2a 2 3a 1 4a 1 5a 1 l7r 3 l^y (2) 

la 2 2a 2 4a 1 5a 2 l^ 2 l7r^ (3) 

where the 5a orbital corresponds to the H 2 lo u orbital. The three configurations 

above were the reference configurations in the POL-CI calculations of Walch et 

•Ms]. 

Table III shows the calculated barrier height and endoergicity obtained using 
CASSCF orbitals for a 3al7r active space and a CCI with eqn. 1-3 as reference 
configurations ([5s4p2dlf/4s2p] basis set). Looking first at the results correlating 
eight electrons, it is seen that the computed barrier height and endoergicity are both 
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too large. Since the saddle point and products regions both contain an 0~ term 
that is not present for reactants, this problem is probably related to the difficulty in 
calculating the electron affinity of O atom at the SCF level. For both saddle point 
and products regions, quadruple excitations, which arise as products of correlation 
terms needed to describe the transfer of the H 2 bond pair and correlation terms 
for describing O - character, are important. These terms may be somewhat less 
important for the products region, since the bond pair has relocalized. As will be 
seen later, this hypothesis is supported by the lowering of the barrier height when 
these higher excitations are explicitly included by expanding the reference space in 
the expanded active space Cl calculations. 

From Table III it is also seen that, at this level of calculation, a better endoergicity 
is obtained if the 02s electrons are not correlated. This is in accord with calculated 
D t ’s for the 0 2 molecule where better results are obtained, when the 02s electrons 
are inactive in the CASSCF calculation, if the 02s levels are not correlated in the 
Cl calculation [22 . As will be seen later, with larger basis sets and a more complete 
correlation treatment, the eight electron endoergicicity is in better agreement with 
the experimental value. For the [5s4p2dlf/4s2p[ basis set, the calculated barrier 
height is the same at the Cl -I- Q level whether six electrons or eight electrons are 
correlated. Thus, the calculations to determine the size of the active space, which 
used this basis set, were carried out with only six electrons correlated. 

Since 0~-like terms are expected to be important in the saddle point region, 
the CASSCF active space was expanded to include the most important of these 
terms. In the present studies the active space was expanded by addition of the 
02p' natural orbitals, where the 2p' is a tight-diffuse correlating orbital (i.e. it 
has an extra radial node like a 3p orbital but is of the same radial extent as a 2p 
orbital), the 03d natural orbitals, and the H 2 2 a g and l7r u natural orbitals, where 
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the orbital notation refers to the reactant geometry. For the a orbitals, the order of 
these natural orbitals (in order of decreasing natural orbital occupation numbers) 
is the same for both the reactants and saddle point geometries: 

H 2 I o g > O2po > H2la u > H22 <?j > 02pa' > 03d 

For the rr orbitals on the other hand the order is different for reactants and saddle 
point geometries. The order for reactants is: 

02p7r > 02px' > H21 tt u > 03d 

while for the saddle point region the order is: 

02p7r > 02p7r' > 03d > H2l7r u 

From the above one sees that in order to include the 02p' shell consistently it is 
necessary to use a ho2T> active space. To consistently add the 03d shell would 
require a 6a47rld active space. 

The effect of adding additional active orbitals was tested at the CCI level by 
systematically expanding the CASSCF active space while adding corresponding 
reference configurations. The smallest active space is 3olir which corresponds to 
the configurations given in eqn. 1-3. The 4 <t27t active space adds the 02p<r' and 
02p7r' orbitals. The 5<r37r active space adds the H 2 2o g like and 03d7r like natural 
orbitals. Finally, the 6cr3rr active space had the same n space as the 5 <t 37T case, but 
the 3dcr like natural orbital was added to the a active space. 

The reference configurations for the MRSD CCI calculations included 1-3 and 
added 4-7. 


la 2 2c 2 3cr 2 4<T 1 27rjl7r* 

( 4 ) 

la 2 2o 2 3cr 1 5<7 1 n<7 1 l7r 2 l7r^ 

( 5 ) 

la 2 2a 2 3a 1 4o 1 5a 1 l7r*/7r*l7r* 

( 6 ) 

lcr 2 2cr 2 3o 1 4a I 5a 1 Inlrmr^ 

( 7 ) 
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In 5-7 n runs from 6 to the total number of a orbitals and / and m run from 2 to 
the total number of 7r orbitals. Eqn. 4 corresponds to the diagonal double 2p— >2p' 
excitation, while 5-7 are interpair terms between the H 2 bond pair and the 2p— >2p' 
excitations. This choice of reference configurations includes all configurations with 
Cl coefficients greater than 0.05 in the CASSCF wavefunction. The barrier height 
was computed in each case with and without Davidson’s correctional] with the 
reactants limit referenced to the SoZt: reactants energy. The energy reference was 
chosen in this way since the reactants energy was relatively insensitive to expansion 
of the active space (the effect of expanding the active space is only 0.3 kcal/mole at 
the reactants geometry) and the resulting values of the barrier height, AEf,, then 
reflect the convergence of the total energy at the saddle point geometry. 

Table IV shows the effect of expanding the active space in the CASSCF calcu- 
lation. Adding the 2p' shell has a large differential effect on the barrier height, 
but addition of further active orbitals has only a relatively small effect. This is 
consistent with the idea that configurations involving the 2p' shell are the most 
important correlation terms needed to describe O - . Based on these calculations 
a 5o2n active space was used in subsequent calculations to determine the saddle 
point properties for O + H 2 . The reference configurations that were used in these 
calculations are given in Table V. The configurations in Table V include a diagonal 
double 3a 2 — ♦ 6<7 2 which is the H 2 ltfp — * > 2 cr 2 and which was not included in the 
set of configurations given above. In addition, the configuration given in eqn(5) for 
n=6 is found to be unimportant and is omitted. 

A grid of points was computed about the expected saddle point and the actual 
saddle point geometry was determined by fitting these data to a function of the 
form: 

F (r j , r 2 ) = c 0 + ci rj + c 2 r 2 + c 3 r 1 r 2 + c A r\ + c h r\ (8) 
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Here the saddle point was determined both for Cl and ,CI 4- Q. At the reactants 
geometry, the H 2 distance was optimized for an 0 4- H 2 supermolecule calculation. 

Table VI shows the results of MRSD Cl calculations for O + H 2 using the 
[5s4p2dlf/4s2p] basis set and a 5a2n active space as described above. The sad- 
dle point geometry obtained here has Ton 0-01 ao shorter and thh 0-06 ao shorter 
than the values determined for the POL-CI calculations. The present studies also 
show a larger curvature in the direction along the reaction coordinate. Based on a 
Wigner correction [23] for tunneling, the larger curvature would indicate tunneling is 
more important on the Cl surface than on the POL-CI surface. The barrier height 
obtained for the 5o27r active space is 2.1 kcal/mole lower than obtained with the 
3ol7r active space( see table III). This is a very significant effect strongly suggesting 
that 0“-like terms are very important for the O 4- H 2 reaction. 

An additional calculation was carried out at the saddle point geometry from Table 
VI using the 5o27r active space CASSCF but with additional reference configura- 
tions. The new reference configurations included: i) selected triple and quadruple 
excitations which arise as products of the diagonal double excitation 3 a 2 — > 5o 2 
and 2p — *• 2p' single and double excitations and ii) atomic interpair terms which 
are double 2p — » 2p' excitations. Inclusion of these addition reference configura- 
tions increases the barrier height very slightly (0.1 kcal/mole). This result indicates 
that the calculation is essentially converged with respect to the list of reference 
configurations. 

Table VII shows two important effects on the barrier height for the O 4 H 2 
reaction as a result of expanding the basis set. The first is that correlating the 02s 
lowers the barrier by only 0.1 kcal/mole with the [5s4p2dlf/4s2p] basis set, but by 
0.5 kcal/mole for the [5s4p3d2flg/4s3pldj basis set. The second effect illustrated in 
Table VII is the importance of higher angular momentum functions for these CPF 
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wave functions. For the case of correlating eight electrons one sees a 1.1 kcal/mole 
lowering of the barrier height for the [5s4p3d2flg/4s4p2dlf] basis set compared to 
the [5s4p2dlf/4s2p] basis set. For the largest basis set the CPF barrier height is 
13.6 kcal/mole when eight electrons are correlated (Table VII). Correcting for the 
basis set superposition error leads to a CPF barrier height of 14.3 kcal/mole. 

Table VIII shows the results of CCI calculations with the ANO basis 6et. Here it 
is seen that the barrier height is 14.7 kcal/mole for CCI and 12.7 kcal/mole for CCI 
+ Q when eight electrons are correlated. For these calculations it was necessary to 
use the contracted Cl, since the calculation involves about 1.2 million configurations 
(uncontracted). For the calculation with the [5s4p2dlf/4s2p] basis set and 5o2ir 
active space the contraction error is found to be 0.5 kcal/mole when six electrons 
are correlated. The superposition error for the ANO basis set is calculated by 
the counterpoise method to be 0.15 kcal/mole when eight electrons are correlated. 
Combining these numbers leads to an estimate of 14.4 kcal/mole and 12.4 kcal/mole 
for Cl and Cl + Q, respectively. Since the multireference Davidson’s correction 
may be an overestimate, we can only say with confidence that the computed barrier 
height is within 2.0 kcal/mole of the currently accepted value of 12.5 kcal/mole; 
although, the true error is probably about half that large. Table VIII also shows 
the computed endoergicity. Here it is seen that the error is 0.9 kcal/mole for CCI 
and 1.6 kcal/mole for CCI + Q. These errors are seen to be of the same order of 
magnitude as the estimated error in the barrier height. 

Table VIII also shows the computed barrier height with respect to the OH + H 
products. The result for CCI + Q correlating eight electrons is 8.2 kcal/mole as 
compared to the experimentally derived value of 9.6 kcal/mole (the 12.5 kcal/mole 
barrier height minus the experimental endoergicity). This result indicates that the 
multireference Davidson’s correction must be overshooting to some extent, since the 
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CCI + Q barrier height is 0.1 kcal/mole and 1.4 kcal/mole lower than ” experiment” 
for the forward and reverse directions, respectively. It should be noted that the 
barrier height for the reverse reaction is not corrected for superposition error or the 
CCI contraction error. 

A remarkable result of these studies is that the present calculations which involve 
more than 1.2 million configurations are finally converging toward the barrier height 
obtained in the POL-CI calculations which involved only a few thousand configu- 
rations. It is clear from the present studies that the POL-CI calculations, which 
allow only one electron outside the 3alx active space, cannot properly describe the 
O' character in the wavefunction at the saddle point geometry, since these studies 
show that a larger active space (5 <t2x) is needed for that purpose. Thus, one might 
have expected that the POL-CI method would result in too large a barrier height as 
was obtained in the CCI studies using a 3ol7r active space. The POL-CI treatment 
neglects the angular and radial correlation of the H 2 molecule which is expected 
be more important for reactants than for the saddle point region. However, the 
neglected ionic terms are less important for the reactants than for the saddle point 
region. Thus, it appears that an accidental cancellation of these opposing effects 
acounts for the good barrier height obtained in the POL-CI. Adding additional 
correlation to this wavefunction by allowing single and double excitations from the 
POL-CI reference set leads to a higher barrier, since angular and radial correlation 
is easy to describe, but a larger active space and extended basis set treatment is 
necesssary to properly describe O" character. It is only after inclusion of extended 
basis sets and active spaces that a balanced description and lower barrier height are 
obtained, as in the present calculation. 

IV. Conclusions. 

The convergence of the barrier height for the O H 2 — * OH -I- H reaction has 
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been studied as a function of the size of the active space and basis set completeness. 
The barrier height is found to be rapidly convergent with respect to expansion of 
the active space. Addition of 2p— >2p' correlation terms to the active space lowers 
the barrier to the O + H 2 reaction by about 2.0 kcal/mole, but addition of 3d 
and other terms has little additional effect. Since the saddle point region of this 
reaction is expected to have a significant 0“ component which will not be important 
for the reactants, these results suggest that the difficulty in accurately computing 
the barrier height for the O + H 2 reaction parallels the problem of describing the 
electron affinity of oxygen. Multireference singles and doubles contracted Cl plus 
Davidson’s correction calculations using a [5s5p3d2flg/4s3p2dlf] basis set with a 
5a27r active space lead to a barrier height of 12.7 kcal/mole. Including an estimate 
of the Cl contraction error and basis set superposition error leads to 12.4 kcal/mole 
as the best estimate of the barrier height. The latter result is in good agreement 
with the current estimate of 12.5 kcal/mole. 
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Table I. Basis Sets for Oxygen. 


O [5s4p] valence basis 


function 

1 

2 

3 

4 

5 

6 

7 

8 
9 

7 

8 

9 

10 
11 


s 

31195.560(0.00021) 
4669.3800(0.00163) 
1062.2600(0.00845) 
301.42600(0.03419) 
98.515300(0.11031) 
35.460900(0.26949) 
13.617900(0.42355) 
5.3862000(0.28304) 
1.5387000(0.02748) 
13.617900(-. 14601) 
5.3862000(-. 14788) 
1.5387000(0.23867) 
0.605500(1.00000) 
0.220500(1.00000) 


P 

64.7719(0.00584) 

14.9727(0.04058) 

4.55440(0.15754) 

1.56370(0.35300) 

0.54107(1.00000) 

0.17776(1.00000) 

0.05840(1.00000) 


[3d2flg polarization basis 


function 

d 

f 

g 

1 

7.5130(0.05799) 

5.4860(0.173786) 

2.5880(0.484830) 

2 

2.3890(0.30456) 

1.9670(0.597338) 

0.9750(0.653938) 

3 

0.9713(1.00000) 

0.8418(1.000000) 


4 

0.4324(1.00000) 




!2dlf; polarization basis 


function 

d 

f 

1 

4.2770(0.16866) 

1.806(0.47694 

2 

1.3410(0.58480) 

0.618(0.65874) 

3 

0.5220(1.00000) 
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Table II. Basis Sets for Hydrogen. 


[4s4p2dlf] Hydrogen basis 

function s pdf 

1 82.63637(0.006172) 6.0760(0.05713) 3.535(0.168660) 1.204(0.476935) 

2 12.40960(0.047210) 1.5750(0.28575) 1.108(0.584798) 0.412(0.658738) 

3 2.823850(0.232530) 0.5554(1.00000) 0.432(1.000000) 

4 0.797670(0.790500) 0.2211(1.00000) 

5 0.258100(1.000000) 0.0880(1.00000) 

6 0.089890(1.000000) 

7 0.030000(1.000000) 


function 

1 

2 

3 


;2pld) polarization basis 


P 

1.7983(0.17705) 

0.4663(0.88560) 

0.1644(1.00000) 


d 

1.827(0.46662 

0.548(0.66447) 
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Table III. Calculations for O + H 2 — ► OH + H 
[5s4p2dlf/4s2p] basis set 
3ol7r active space 


CCI a,fc , kcal/mole 


Geometry 

6 electrons 

8 electrons 

expt. 

products 

3.6(3. 3) 

6. 7(6.6) 

2.9 

saddle point c 

17.0(16.0) 

17.5(16.0) 

12.5 

reactants 

0.0 

0.0 

0.0 

De(H 2 ) 

4.62 

4.62 

4.75 

De(OH) 

4.48 

4.33 

4.62 

0 Values in parenthesis include Davidson’s correction. 
b All quantities are in kcal/mole except the D e values which 

are in eV. 

c The saddle point geometry is tqh 

=2.293 and r H H = 1.764 

which is very close to 

the POL-CI saddle point geometry. 

The OH and H 2 bond lengths are tqh = 1.868 


and Thh = 1-431 
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Table IV. Calculations for O + H 2 — ► OH + H 
[5s4p2dlf/4s2p] basis set 
Effect of expanded active space 


AEb °, kcal/mole 


active space 

CCI 

CCI + 

3crl7r 

17.3 

15.9 

4a27r 

15.2 

14.4 

5<t37t 

15.1 

14.3 

fkr37r 

15.0 

14.2 


0 Using the reactants and saddle point geometry as in Table III. 
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Table V. Reference Configurations for O + H2. 


l<r 2 2cr 2 3cr 2 4o 1 l7r^l7r^ 
lo 2 2a 2 Za 1 4a 1 5o 1 l7r 2 l7r* 

la 2 2fT 2 4o 1 5a 2 l7r 2 l7Ty 
1 o 2 2 o 2 4 a 1 60 2 1 7T 2 1 7T * 
l(T 2 2cr 2 3o 2 4a 1 27r 2 l7r^ 
la 2 2<T 2 3£r 1 5a 1 7a 1 l7r 2 l7rJ 
lcr 2 2a 2 3o 1 4a 1 5<7 1 l7r^27r]l7r^ 
la 2 2<T 2 3o 1 4a 1 5a 1 l7r 2 27r’ 
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Table VI. Computed Saddle Point Properties for 0 + H 2 . 0 


property 

POL-CI 

Cl 

CI + Q 

AEb, kcal/mole 

12.5 

14.9 

13.9 

r on, ao 

2.320 

2.313 

2.314 

ao 

1.747 

1.690 

1.688 

curvature b 

-0.057 

-0.083 

-0.078 


0 Multireference SDCI calculations with a 5a2?r active space and a [5s4p2dlf/4s2p] 
basis set. 

** Negative eigenvalue of the Hessian matrix. 



Table VII. Calculations for O + H 2 - OH + H 
Effect of Basis Set 


AEb °, kcal/mole 


Basis set 6 electrons 8 electrons 

[5s4p2dlf/4s2p] 14.83 14.72 

[5s4p3d2flg/4s3pldj 14.73 14.20 

[5s4p3d2flg/4s4pldj 14.09 

[5s4p3d2flg/4s4p2d; 13.87 

[5s4p3d2flg/4s4p2dlf] 13.60(14.3 b ) 


0 CPF calculations using the saddle point and reactants geometries from the mul- 
tireference SDCI calculations using a 5o2n active space and the [5s4p2dlf/4s2p] 
basis set. 

Including an estimate of basis set superposition error. 
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Table VIII. Computed Barrier Height and Endoergicity for O + H 2. 0 


property 
AEj, b , kcal/mole 
contraction error 
superposition error 
AEt,(corr.) 

AE rs , kcal/mole 


CCI(6 electrons) 
15.1(14.0) 

0.5 

( 0 . 0 ) 

14.8(13.7) 

1.7(1. 9) 


CCI(8 electrons) expt. 

14.7(12.7) 

0.5 

0.15 

14.4(12.4) (12.5) 

3. 8(4. 5) 2.9 


0 Multireference CCI calculations with a 5o2ir active space and a [5s5p3d2flg/4s3p2dlf] 
A NO basis set. 

fc The computed barrier height viewed from the OH -f H direction is 13.5(12.1) and 
11.0(8.2) kcal/mole for correlating six electrons and eight electrons, respectively. 
This may be compared to 9.6 kcal/mole for * experiment”. 
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